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Two targeting schemes have been known for the density matrix renormalization group (DMRG) 
applied to non-Hermitian problems; one uses an asymmetric density matrix and the other uses 
symmetric density matrix. We compare the numerical efficiency of these two targeting schemes 
when they are used for the finite temperature DMRG. 
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I. INTRODUCTION 

The density matrix renormalization group (DMRG) es- 
tablished by White |jj has been successfully applied to 
various problems in condensed matter physics. Q A re- 
cent technical progress in DMRG is its applications to 
non-Hermitian problems, such as asymmetric exclusion 
process, ||| I reaction-diffusion process, || and quantum 
Hall effect. |§ 

For these non-Hermitian problems, left eigenvectors of 
the Hamiltonian are not equal to the complex conjugates 
of the right eigenvectors. Two different targeting schemes 
have been used for DMRG under the situation. One is 
to use an asymmetric density matrix, which is a par- 
tial trace between the left and the right eigenvectors. |7]] 
This scheme has been used for the DMRG applied to 
classical systems || and the finite temperature (finite-T) 
DMRG. |J |ll| The other scheme is to use a symmetric 
density matrix, which is created by targeting both left 
and right eigenvectors as two individual vectors. |Q,||,||j6) 

The purpose of this paper is to compare the numerical 
efficiencies of these two schemes, by observing the cut-off 
error of the renormalization group (RG) transformation 
applied to the finite temperature Heisenberg spin chain. 
In the next section we define the cut-off errors as a func- 
tion of projection operator, that represents the freedom 
restriction by the RG transformation. The two target- 
ing schemes are briefly reviewed in §3, and these schemes 
are compared by calculating the cut-off error numerically. 
Conclusions are summarized in §4. 



II. CUT-OFF ERROR IN RG 
TRANSFORMATION 



system; Z = TrT N . When N is sufficiently large, Z is 
well approximated as 
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where V L and V R is, respectively, the left and the right 
eigenvector of T, that satisfies the eigenvalue relation 
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We have used the normalization (V L , V R ) = £V • V^V R 
= 1 in Eq.l. It should be noted that V L is not equal to 
V R in general, because of the asymmetry T ^ T T . 

Let us consider a formal decomposition of these vectors 
into the products of matrices 
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according to the convention in DMRG, [0 where L / R 
and Q L / R satisfy the orthogonal (or duality) relations 
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The finite-T DMRG estimates the free energy 

of one-dimensional quantum systems by way of a precise 
approximation for the largest eigenvalue A of the quan- 
tum transfer matrix |Q (QTM) T. As an example of 
QTM, we consider that of the 5=1/2 Heisenberg chain. 
(See Fig.l.) We express the matrix element of the QTM 
as T V y ij , where i {%') and j (j 1 ) represents the upper-part 
(U-part) and the lower-part (D-part jL3|) of the column 
spin, respectively. It has been known that the partition 
function of TV-site system can be approximated by that 
of the Trotter |Q decomposed two-dimensional classical 



and £ and rj represent the block spin variables. The 
matrices L / R and Q L / R play the role of renormaliza- 
tion group (RG) transformations, when the freedom re- 
striction 1 < £, 77 < to is considered for for both V L 
and V R in Eq.3. Under the restriction, Vp and vfi are 
m-dimensional matrices that represent the renormalized 
states. 

In DMRG applied to classical system or finite tempera- 
ture quantum system, the RG transformations L / R and 
ql/r are determined so that the cut-off error of the par- 
tition function 
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FIG. 1. Quantum transfer matrix T^, i - of the 5=1/2 
Heisenberg spin chain, where the Trotter number M is 3, and 
both i (i') and j (j') consists of three S = 1/2 spin variables. 

5Z = Tr (1 - P)T N = (V L , (1 - P)V R ) \ N (5) 

is suppressed, where P is the projection operator 

p _ pU pD 

that represents the Hilbert space restriction by the RG 
transformation; P u and P D is projection operator for U- 
and D-part, respectively. Since the operator 1 — P = 
1 — p u p D in Eq.5 can be factorized as 

(1-P U ) + (1-P D )-(1-P U )(1-P D ), (7) 

and the third term is negligible when it is applied to V L 
and V R , we can precisely estimate the relative cut-off 
error 6Z/Z by calculating the inner product 

(V L , (1 - P U )V R ) + (V L , (1 - P D )V R ) 

= (1-TrPV) + (l-TrP D j0 D ) , (8) 

where p u and p D are the asymmetric density matrices 

3 

!>}., E^* (9) 

i 

that satisfy the normalization Trp u = Trp D = (V L , V R ) 
= 1. Let us keep in mind that TrP u p u and TrP D p D are 
essential for the cut-off error of the RG transformation 
by P = P U P D . 



III. ASYMMETRIC AND SYMMETRIC 
TARGETING 

Two different targeting scheme have been used to de- 
termine the RG transformation matrices O l / r and <5 L ^ R . 



One is to obtain them by diagonalizing the asymmetric 
density matrices @,|,0[n) in Eq.9 

( 10 ) 

where is the common eigenvalue for both p u and p D 
in the order of decreasing absolute value; normally all the 
are positive. The projection operators P u , P D , and 
P created from L / R and Q L / R in Eq. 10 are asymmetric, 
as p u and p° are. Let us call such a construction of P u / D 
as symmetric targeting. In this case, the relative cut-off 
error in Eq.8 can be calculated from the eigenvalues of 
the asymmetric density matrix as 

m 

2-TrP u p u -TrP D p D = 2(1-^^). (11) 

? 

It is possible to choose 0^ R and Q\^ K so that vf? and 
v R v become simultaneously diagonal: vjL — v R v = 5^ oj^ 
where w| — . Thus we can interpret the decomposition 
in Eq.3 as an extension of the singular value decomposi- 
tion for the dual vectors V L and V R . Q 

The other targeting scheme is to treat V L and V R 
as individual vectors, |^|| as they simultaneously target 
ground and excited states in DMRG applied to Hermi- 
tian quantum systems. B In this case, the RG trans- 
formations are obtained by first creating the symmetric 
density matrices 

3 3 
i i 

and then by diagonalizing them 
n 

where we have dropped the label L and R from O l / r and 
Q L/R , because O h = O n and Q L = Q R . In this case, the 
projection operators 

m 

Pi'i = £ Oi'tPit 
m 

P?3=Y.QrnQ 3 n (W) 
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are symmetric. Let us call such a construction of P u / D as 
symmetric targeting. Unlike the asymmetric targetting, 
it is impossible to make the m-dimensional matrices 
and simultaneously diagonal. This targeting scheme 
is often used because there is no need to diagonalize the 
asymmetric density matrix, which requires special nu- 
merical care. |l5) It should be noted that the relative 
cut-off error 
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in the symmetric targeting is not directly related to the 
eigenvalues of the symmetric density matrices in Eq. 12. 

Now let us compare the relative cut-off errors for both 
asymmetric and symmetric targeting. We choose the 
5 = 1/2 isotropic Heisenberg spin chain as the refer- 
ence system, and fix the Trotter number M = 7 so that 
we can obtain the eigenvector of T; we have to know V L 
and V R exactly in order to evaluate 6Z/Z. We consider 
the case where U-part contains the same number of spin 
variables as D-part; this U-D division is normally used 
for the infinite system algorithm, [jjj Since U-part is iden- 
tical to the D-part, 1 - TrP u p u _= 1 - TrP D p D holds for 
Eq.ll and l-TrP u p u = l-TrP D p D for Eq.15. Figure 2 
shows the relative cut-off errors when the imaginary time 
step J At O is equal to 1/7. As it is seen, 1 - TrP u p u 
decreases monotonically with respect to m, and is always 
positive. On the other hand, the dumping of 1 — TrP u /? u 
with respect to m is oscillatory; 1 — Tr P u p u is not al- 
ways positive, ||l6|| and the calculated partition function 
is not the variational lower bound. For most of m the er- 
ror 1 - Tr P u p u is smaller than 1 1 - Tr P u p u | , that shows 
the superiority of asymmetric targeting for the finite-T 
DMRG. 
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FIG. 2. The relative cut-off error for the asymmetric tar- 
geting 1 — TrP u p u in Eq.ll and that for the symmetric tar- 
geting 1 - TrP u p u in Eq.15 when J At = 1/7. For the lat- 
ter we use triangle mark when it is positive, and use cross 
mark when negative. The horizontal axis shows the number 
of states kept. 
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FIG. 3. The relative cut-off errors in relatively high tem- 
perature J At — 1/14. 

Figure 3 shows the cut-off errors for relatively high 
temperature J At = 1/14. In both Figs. 2 and 3, we have 
to keep twice as large as m for the symmetric targeting 
in order to keep the same cut-off error of the asymmetric 
targeting. This may be explained by the fact that the 
asymmetric projection operator P u / D is created by 2m 
numbers of linearly independent vectors, while the sym- 
metric projection operator P u / D is created by m numbers 
of orthogonal vectors. 



IV. CONCLUSION 

We have compared the numerical efficiency of the sym- 
metric and asymmetric targeting schemes when they are 
applied to the finite temperature DMRG. It is shown that 
the cut-off error calculated by the symmetric targeting is 
larger than that of asymmetric targeting; as far as cut-off 
error is concerned, the asymmetric targetting is superior 
to the symmetric targetting. 

If we keep twice as large as m for the symmetric target- 
ing, we can recover the numerical precision of the asym- 
metric targeting. Therefore, for the problems that does 
not require large m for the DMRG calculations, the sym- 
metric targeting is of use, in the sense that it does not 
require the diagonalization of asymmetric density matrix, 
and is free from complex eigenvalue problem. [0,fl5| 
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